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H ' Abstract 
■ 

The \^D2 (cc) state is the ground state of spin-smglet D-wave charmonia. Al- 
\^ . though it has not been found yet, the experimental data accumulate rapidly. This 

charmonium attracts more and more attention, especially when the BaBar Collabo- 
(— I ■ ration finds that the A (3872) particle has negative parity. In this paper we calculate 

0^. the double- gamma and double-gluon annihilation processes of charmonia and 

p) bottomonia by using the instantaneous Bethe-Salpeter method. We find the rela- 

O ■ tivistic corrections make the decay widths of \^D2 (cc) 2~5 times smaller than the 

non-relativistic results. If this state is below the D^D^ threshold, we can use the 
sum of annihilation widths and EM transition widths to estimate the total decay 
width. Our result for (cc) with m = 3820 GeV is T = 432 keV. The dominant 

on ' decay channel 1^D2 (cc) — t- /ic7, whose branching ratio is about 90%, can be used 

, to discover this state. 

O ■ 1 Introduction 
m 



> 
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The 2~+ state draws more and more attention [1, 2, 3, 4, 5, 6, 7, 8] recently. One main 
reason is that the latest result of the BaBar Collaboration [9] favors negative parity 
for A(3872), which has induced many discussions about the possibility of assigning the 
1^-D2 (cc) state to this particle. However, theoretical calculations show that the 2""*" 
assignment strongly contradicts with the experimental data of the electromagnetic decays 
of A(3872) [2, 3, 4, 8]. It's very interesting to study the properties of the 1^D2 (cc) state, 
which is the only ground state of spin-singlet D-wave charmonia. The discovery of this 
particle will greatly support the quark potential models and be helpful in understanding 
the non-perturbative properties of QCD. 

These models predict the mass region of this particle is 3760 ~ 3840 MeV [2], below 
the threshold value of D^D^*. This means no OZI decay channels are allowed. So the 
electromagnetic and light hadronic decay channels are important for the discovery of 
this particle. For example, the clean diphoton decay channel will play an important 
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role in determining the inner structure of these particles. This channel can be used 
to distinguish mesons with a qq structure from those without [12]. We note that non- 
relativistic calculations of the two-gluon annihilation process of D-wave mesons have been 
performed recently in Refs. [1, 10, 11]. All of them get a relatively large result. In 
Ref. [1], Chao et ai. have calculated this process to the third order of within NRQCD 
formalism. They show the next-to-leading order QCD corrections contribute enhancement 
factor of 1.8 and 1.5 for charmonium and bottomonium 1^D2 states, respectively. With 
the non-relativistic approximation, the decay widths of D-wave mesons are related to the 
second derivative of the radial wave functions at the origin. This method will result in 
large errors, since the full behavior of the wave function (or the relativistic correction) is 
significant for D-wave mesons. 

The semi-relativistic calculation [12] and the relativistic calculation [13] have been 
performed before. In Ref. [12], the transition amplitude is calculated based on the non- 
relativistic wave function and free quark-antiquark annihilation Feynman diagrams, and 
the dependence of the meson mass is introduced by adding an additional term in the 
Lagrangian. Ref. [13] gives a relativistic result, but there they use a potential with 
the timelike vector spin structure and expend the wave function (amplitude) in a set of 
Laguerre basis functions. Since large relativistic correction is expected for D-wave meson, 
a careful relativistic calculation with the wave function given by a more reasonable way 
is very necessary. This will be helpful for the discovery and study of this particle for the 
future experiments. 

In the previous papers [14, 15, 16], we have calculated the two-photon and two-gluon 
annihilation rates of '", 0"^"'" and 2^^^ states by using the relativistic Salpeter method [17, 
18]. Large relativistic corrections are found, especially for the two P-wave states. There we 
made an approximation that the time components of the quark-antiquark momenta inside 
a meson are constants, setting = P2 ~ T (R-^fs. [12] use the similar approximation). In 
this paper, we study the annihilation processes of 2 ^ states within the same formalism, 
but with some improvements. Except for the method we used before, we also perform the 
calculation without that approximation and compare the results of these two methods. 
We examine the accuracy of this approximation and provide some useful information for 
future study. 

The paper is organized as follows: In section 2, we present the general form of the 
2^+ wave function and the coupling Bethe-Salpeter (BS) equations fulfilled by this wave 
function. We also give the transition amplitude for the two-photon (gluon) annihilation 
processes within Mandelstam formalism. In section 3 we show and discuss the numerical 
results. 



2 Theoretical calculations 

The general wave function of the 2 ^ state with mass M , momentum P and polariztion 
tensor e^i, can be written as [8] 
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Figure 1: Feynman diagrams for 2 — 77. 

where q'^ = — j^P^; q is the relative momentum between the constituent quark and 
antiquark. The constraint condition (see the last one in Eq. (13)) of the scalar functions 
fi has the following form: 



/3(g±) 



fi{q±)M{miUJ2 - m2UJi] 

q'iiui + U2) 
-f2{q±)M{ui + U2) 



(2) 



{miU2 + coim2) 

where cuj = rrif — q± = \/m'^ + (p. In the equal mass case, we get cui = UJ2 and /s = 0. 

With the same method used in Ref. [19], we can get the coupled instantaneous BS 
equations for the 2~~^ state 

(M-2u;i)(A(g-) + ^/2(g1) 
mi 

d'k—^ [(g ■ kf - \(fk'][m,iVs - V,){m,f2{k) + ujrh{k)) - (K + V,)k ■ qf2{k)l 

Zq TTliUli o 

{M + 2uj,)UM-—f2m 
mi 

d'k—^ [(g ■ kf - ^(fk'][rm{Vs - V,){mif2{k) - Uih{k)) - (K + V,)k ■ qf2{k)]. 

(3) 

By solving above equations, we can get the mass spectrum and corresponding wave 
functions, /i and /2 fulfill the normalization condition: 



d?q 16 ^ ^ Wi _4 



(27r)3 3 ''•'mi 



hf2—<t = lOM. 



(4) 



According to Mandelstam formalism [20], the relativistic transition amplitude for the 
double-photon decay processes (see Fig. 1) can be written as 



(mi — ie) 



'1 ' ri, 



\-h- i^i - ^e)' 



Mi+Mr 



J} 



(5) 



where xil) is the BS wave function of the initial meson; ki and are momenta and 
polarization vectors of final photons, respectively; = | for charmonium and = — ^ 
for bottomonium. 
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Figure 2: Wave functions of l^D2{cc) and rjc {I^Sq). 



In the following, we give the detailed calculation for Aii in Eq. (5). 



(mi — ie) 



\/3(ee,)^ 



y/3{eeq 



Tr{x(g)- 



Tr|x(g) 



aiP + qY - 2(aiP + g) ■ fci - mf + ie 



} (6) 



(27r 



■}• 



In the second line we have used p^* = aiP^ + g'^, where the third line we 

have used = q^^ + q^ and q^^ = jjP'^, where Qp = We also define 



= -{ql - — -2q^-h- ml) = {q- kif + m^. 



(7) 



The BS wave function xio)-, which is related to Eq. (2) through 
has the form [19]: 



where 



X(g) = Si{pi)T]{qA_)S2{p2), 

v{q±)= d^k±V{q±,k±)ip{k±). 



(9) 
(10) 
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Si{pi) is defined as 



qp + aiM — wi + ie qp + aiM + Ui — ie' 



(11) 



-^2(^2) = ■ 77 — + 



where 



—qp + — UJ2 + ie —qp + + 002 — ie' 



At = ^^[-u.2Ti^, + m2)]. 



(12) 



When mi = 772,2, we have wi = U2 and A]*^ = A2'. 

Inserting Eq. (9) and Eq. (11) into Eq. (6) and performing the counter integral over 
qp by considering [19] 

{M-uji- uj2)^+-^{q±) = A^(g±)r/(g±)A^(g±), 

(M + c^i + uj2)v~-{q^) = -A7(gx)r/(g±)A^(gx), (13) 
f^-{q±) = V'^{q±) = 0, 



we get 



where 



M, =zV3{ee,y J ^{^iTrb++^((l - + i^-h + rn,)i,] 

+ C2Tr[v,++^2((ai - + i^-h+ ^1) A] 

+ C3Tr[^-/,(-^/' + - ^1 + mi)/J 

+ C4Tr[(/.-/2((«i - + i^-h+ ^1) A] 

+ C5Tr[A+7/(g^)A2-/2((«i " + i^~h + ^O/J 

+ C6Tr[Arr/(gx)A+^2((«i + + " ^1 + "^1) All' 



_ 1 ^ M-2a;i 

'-'1 ~ ~7 ;r7 v5 To' '-^2 



(14) 



(aiM-u;i)2-A2' " 2A[(aiM-wi)2- 

1 ^ _ M + 2a;i 

' (aiM + a;i)2-A2' ^ 2A[(aiM + a;i)2 - ^2] ' ^ ^ 

^5 = ^*6^ 



2A[(A + a;i)2- (aiM)2]' 

ip^^ and </9 are the positive and negative parts of the wave function, respectively. In 
the following, we will just consider the contribution of ip^^ . This term is dominant and 
other terms can all be ignored. 

ip^^ has the following form: 

<^++ = e^.q1.ql{A, + + A^f(^J^^, (16) 
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where 

A = ^(/l + ^/2), ^ = ^(/l + — /2), ^3 = -^(/l + — /2). 

2 mi ZMuji mi IMuji mi 

Inserting Eq. (16) into Eq. (14) and finishing the traces, we get 

Ml = -AV3iee,)ha(Sase'^e^,P'k'[e^, J ^qIqUCi + Cs)^. 

The integral can be expressed as the combination of g'^'^, ki and P. Through 
calculations, we get 

Ml = V?>{ee,feo.p.5e1e^2P'^^Xk'iKE. 

E is defined as 



where 9 is the angle between q and ki. 

With the same method, M2 can be expressed as 

M2 = V^{eeqfeo.pase1e^^P'e^,k'tk\klE\ 

where the prime means that we have changed /ci, ei into /c2; ^2 in Eq- (6). 
By adding Eq. (19) and (21), we finally get the transition amplitude 

M = Mi + M2 = y3(ee,)2e„^.5e?efPV^r^i^r + E'). 
With the completeness condition of the polarization tensor 

A 

p p 1 p p p p 

we get the unpolarized transition amplitude squared 

The two-photon and two-gluon annihilation rates have the respective forms: 
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3 Results and discussions 



When solving the BS equation, we have used the Cornell potential 

^(g1 = K(g)+7o®7°K(g-), 



A 



a 



IT- 



a 



2\2 ■ 



127r 



2 as{q) 
1 



(26) 



27 \n{c 



^QCD ' 



2.7183, 



a 



0.06 GeV, 



where the following values for the parameters are used: a = ( 
X = 0.21 GeV^ = 1.62 GeV, mb = 4.96 GeV, Aqcd = 0.27 GeV (for bb, Aqc7d=0.20 
GeV). Vq can be fixed by fitting the mass spectrum. However, there is no experimental 
value available now. So we just take M{lD{cc)) = 3820 MeV as an example which 
lies in the range of 3.760 ~ 3.840 GeV predicted by quark potential models [2]. As for 
M{lD{bb)), we take 10.15 GeV which is the same as that in Ref. [21]. By doing so, we 
get Vq = -0.144 GeV for charmonium and -0.15 GeV for bottomonium. 

We present masses of the first five states in Table 1. One notices that for charmonium 
only the ground state is below the DD* threshold, while for bottomonium the first two 
states are below the BB* threshold. For these states, we can use the sum of double-gluon, 
double-photon and El decay widths to estimate the total width. 



Table 1: The predicted masses (GeV) of 2 + charmonia and bottomonia. 



Meson 


ID 


2D 


3D 


4D 


5D 


cc 


3.820 


4.151 


4.405 


4.611 


4.781 


bb 


10.15 


10.45 


10.70 


10.90 


11.08 



Table 2: Double-gamma decay widths (eV) of 2 charmonia and bottomonia. Results 
in parentheses are obtained by using the Pi=P2 = y approximation. 



Meson 


ID 


2D 


3D 


4D 


5D 


cc 


11.6 (14.8) 


13.4 (18.7) 


13.4 (19.4) 


12.9 (18.9) 


12.0 (17.9) 


bb 


0.0475 (0.0590) 


0.0768 (0.0959) 


0.0958 (0.120) 


0.108 (0.135) 


0.116 (0.146) 



Except for the above method, we also perform the calculation as we did in Ref [14, 15, 
16]. That is, the dependence of the integrand in Eq. (5) only comes from x(g) with the 
condition = ^ (which is equivalent to pi = aiP + q±). Using Eq. (8), we can integrate 
out g° and get 



M = -i 



{2n) 



1 



1 - ^1 

d^q± . 



[mi — le 
1 



1 



{2Tif\q + kiY 



+ 



mi — le) 
)qlqlf2. 



(27) 
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Table 3: Double-gluon decay widths (keV) of 2 charmonia and bottomonia. Results in 
parentheses are obtained by using the = P2 = T approximation. 



Meson 


ID 


2D 


3D 


4D 


5D 


cc 


35.6 (45.6) 


41.3 (57.5) 


41.4 (59.6) 


39.6 (58.1) 


37.1 (55.1) 


bb 


0.883 (1.20) 


1.43 (1.78) 


1.78 (2.23) 


2.01 (2.52) 


2.16 (2.71) 



Table 4: Double-gamma decay widths (eV) of 2 charmonia and bottomonia of different 
models. In parentheses, meson masses (GeV) of different models are presented. 



Model 


Ours 


Ref [12] 


Ref [13] 


Ref [23] 


Ref [22] 


ID(cc) 


11.6 (3.820) 


20 (3.84) 


13.6 (3.84) 


65.7-73.2 (3.872) 


~ 50 (3.5 ~ 3.7) 


2D{cc) 


13.4 (4.151) 


35 (4.21) 


20.2 (4.28) 






lD{bb) 


0.0475 (10.15) 


0.033 (10.15) 


0.0513 (10.13) 






2D{bb) 


0.0768 (10.45) 


0.069 (10.45) 


0.0962 (10.47) 







Table 5: Double-gluon decay widths (keV) of 2 charmonia and bottomonia of different 
models. In parentheses, meson masses (GeV) of different models are presented. 



Model 


Ours 


Ref [1]* 


Ref [10] 


Ref [24] 


Ref [11] 


Ref [22] 


ID(cc) 


35.6 


155 


110 (3.815) 


60 


190 (3.837) 


~60 


IDibb) 


0.883 


3.22 











* This is the leading order results. To the 0{al) order, the results is 274 keV and 4.70 keV for cc and 
hh, respectively 



One notices that this result has a similar structure as in Eq. (18). Actually, if we make 
the approximation mi ~ wi ~ aiM and /i ~ /2, Eq. (18) will reduce to Eq. (27). 

In Table 2 and 3, we give the two-gamma and two-gluon decay widths for the first five 
2 ^ heavy quarkonia. The results of two methods are in good agreement with each other, 
while the one with Pi = Y approximation is larger. As the principal quantum number 
increases, the decay width first increases and then decreases for charmonium, while it 
increases all the time for bottomonium. This can be understood as follows. First, the 
overlap integral E {E') in Eq. (20) decreases as the principal quantum number increases 
(for bottomonium, E first increases, then decreases). This is because the cancellation is 
more severe if the wave function has more nodes. Second, there is also a factor M which 
is increasing all the time. The final result comes from combining these two effects. 

We present two-gamma decay widths for the first two charmonium and bottomonium 
states with different models in Table 4. Our result is close to that of Ref. [12] and [13], 
which are based on the semi-relativistic and relativistic formalism, respectively. Ref. [22] 
gives the non-relativistic result which is about 4 times of ours. With the two-body Dirac 
equation (TBDE) method, Ref. [23] gets a much larger result. In Table 5, two-gluon decay 
widths for the ground state with different models are given. Ref. [24] use the known values 
of charmonia decay widths into different channels and the non-relativistic formulas to get 
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a result which is about 1.7 times of ours. Ref. [10, 11, 22] uses the non-relativistic method. 
The deviation of the results comes from the use of different parameters. Using the NRQCD 
method, Ref. [1] gets the leading order result which is about 4 times of ours. 

Our result shows great relativistic corrections exist in the annihilation processes of 
D-wave quarkonia. This can be understood from the behavior of the wave function in 
momentum space. In Fig. 2, we plot the ground-state wave functions of the and 
charmonia. For l^D2{cc), one can see the wave function reaches its maximum when q 
takes a relatively large value (~ 0.7 GeV), which means the relativistic corrections will 
be significant. This can be compared with the rjc (S-wave) case. There the dominant 
contribution of the wave function comes from the range where q is small, but even so the 
relativistic corrections is about 30 ~ 40%. 

For the state below DD* or BB* threshold, we can estimate its total decay width by 
adding the widths of different channels together as mentioned previously. For l^Z)2(cc), 
we can get F ^ Vgg + F^^^^ + F^j/^ = 432 keV, where we have used the results of F^^^ and 
F^j/^ in Ref. [8]. With this total decay width we can estimate the branching ratio of the 
two-gamma channel: Br(l^D2 (cc) — )■ 77) = 2.7 x 10~^. This is close to that of ?7c, while 
is one order of magnitude less than that of Xco and Xc2 [25]. The golden channel to find 
the 1^D2 (cc) is its decay to 7/ic since its branching ratio is about 90%. 

In summary, by using the instantaneous BS method, we have calculated two-photon 
and two-gluon annihilation processes of charmonium and bottomonium states. It shows 
that the decay width of the charmonium increases first and then decreases, while the 
decay width keeps increasing for the bottomonium case. The relativistic corrections are 
very large, which implies one has to include them. Although the two methods we used 
in this paper give results very close to each other, the one without the = P2 ~ T 
approximation is more justifiable. 
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